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A. Objectives 

The objective of this research is twofold: 

1. Testing of an equilibrium hypothesis for the density of 
natural vegetal systems. This is being accomplished by comparing data from 
a wide range of climate and soil conditions to a hypothesized physical rela- 
tionship between the vegetal canopy density, and the ratio of actual to 
potenti.al evapotranspiration. 

2. Use of this hypothesis in an existing mean annual water balance 
equation to determine the spatial average effective soil properties for 

a given catchment. For verification, these soil properties will be used 
in the water balance equation to derive the statistics of annual watershed 
yield which will be compared with the observed statistics. 

B . Technics ‘ ogress 

1. Testing of Hypothesis 

Based upon survival-of-the-f Ittest arguments, Eaglcson [1978] has 
developed a theoretical relationship between the ratio, J, of actual to 


potential evupotransplration (l.e., the evaporation efficiency) and the 

vegetation canopy density, under natural cquilibrlua conditions. Ttiis 

relationship is presented in Fig. 1 for varlci.s values of the parameter 

6h /e . Here 
o p 

6 ^ ■ expected value of the t'.me between rainstorms, days 

h^ > depth of water retained on the surface, cm 

e • average potential rate of evapotranspiration from a 
P 

bare soil surface 

Also shown on this figure are pr ats representing a preliminary 
analysi*' of observations from catchments in various climatic regimes 
throughout the U.S. 

The two points labeled Santa Paula and Clinton are the two most 
thoroughly studied watersheds. Their equilibrium vegetal densities have 
been estimated at .3 and .8, respectively. Because of the in-de-pth analysis 
of the climate and watershed conditions of these two catchments, all the 
values of necessary parameters and variables were established with little 
uncertainty, and it can be se^'n that they afford the best fit to the derived 
rel. 1 t ionshlp. 

All points labeled "W - ***", are catchments whose precipitation 
and runoff data and general watershed conditions were studied and published 
by the U.S.D.A. [1963, 1967j. Difficulties arose with the sets of data in 
that many of the catchments were not in their natural state thus resulting 
in vegetation densities that may not be the natural equilibrium density. 
Also, estimates of the albedo and cloud cover, which are important in the 
calculation of potential evaporation, are not exact. The two points labeled 
W-3 and W-2 are catchments in Florida where groundwater infiltration from 


neighboring wacershede Is highly possible due to the flat nature of the 
topography. This infiltration could increase the actual yield of the basin, 
thus decreasing the observed evapotransplration, resulting in the point 
being located below the derived curve. Efforts are being made to obtain 
more exact data on these watersheds. 

The points marked Boco Mtn. and Palo Al'.o were obtained from a 
paper by F. A. Branson and J. B. Owen (19781. In this report, the authors 
plotted a relationship between runoff and percent bare soli in a watershed. 
Also given In the paper was the mean annual precipitation for these two 
catchments. Potential evaporation was calculated for these catchments, 
for those studied by U.S.D.A. , and Tor Clinton and Santa Paula, from mean 
annual values for temperature, cloud cover and relative humidity obtained 
from U.S. Weather Service publicatloas for the nearest weather station. 

The evapotransplration efficiency, J, was estimated by subtracting 
annual stieamflow from annual precipitation, and dividing by the potential 
evaporation, calculated as explained above. 

2. Estimation of Effective Soil Parameters 

To estimate the spatial average effective soil parameters, K(l), 
the saturated intrinsic permeability; c, pore disconnectedness index; and 
s^, average soil moisture concentration, the following algorithm is being 
used: 

a. For a given catchment, the vegetation parameters, and K^, 
are known from observations 

b. A value for n, effective soil porosity, is assumed 

c. From the relaticnship between M^, J and E, and using the 



3J 

the hypothesis that * 0 at equilibrium, the value of E 

can be calculated from M - M . This gives Equation (1) 

o 

S - 8 (K(l), c) (1) 

o o 

d. Extrapolating the cuive of J vs. E from the knunai value of £ 
at the given to the point J * 1, by some sort of 3rd order 
approximation, the value of E(l) which corresponds to s^ ■ 1 
can be determined from the function [see Figure 2] derived by 
Eagieson [1978]. This gives Equation (2) 


s - 1 - 8 (K(l), c) 

O 3 


( 2 ) 


e. The equation for E is 
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where K(l), <|>(1), d, are functions of K(l) and c, and 

the other variaoles are known, for s ■ 1, E ■ E(l), /.nd the 

o 

above equation can be aclved directly for K(l) with c as the 
only independent variable. 

f. It can al^o be seen from Equation (3) chat the only difference 
between E(l), and E is the last term. With E(l) and E known, 
and for n given, c and can be directly calculated from 


s 

o 


iKd) 


]/d+2 


( 4 ) 


g. Picking a value for c, K(l) and can be calculated. The 
last eq>.:tion needed is the water balance equation. 





( 5 ) 


8 - 8 (K(l), C) 

o o 

A trial can Chen be performed on c until the water balance 
equation ia aatiafied. 

Preaently, difficulties are being experienced in the approximation 
for E(l). Initially, a simple tangent approximation to Che point J 1 
from the known value of E was used to find E(l). However, due to Che 
asymptotic behavior of the function in Che area around J - 1, this approach 
resulted in values of E(l) which are smaller than the actual values. These 
actual values were calculated from data obtained at Clinton and Santa Paula, 
and were found to be two orders of magnitude greater chan those obtained 
from the tangent approximation. Use is now being icade of a Taylor series 
expansion of higher orders to simulate the asynptotic behavior of the func- 
tion. 
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